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Abstract. We present results of an extensive investigation of the physical changes taking place in the emission 
regions of the LMXB X 1624-490 during strong flaring in observations made using RXTE in 1997 and 1999. Based 
on the detailed light curve, we propose that the flaring consists of a superposition of X-ray bursts. It is shown 
that major changes take place in the blackbody emission component, the temperature kTsB increasing to ~2.2 
keV in flaring. Remarkably, the blackbody area decreases by a factor of ~5 in flaring. During flare evolution, 
the blackbody luminosity remains approximately constant, constituting a previously unknown Eddington limiting 
effect which we propose is due to radiation pressure of the blackbody as /cTbb increases affecting the inner disk 
or accretion flow resulting in a decreased emitting area on the star. We argue that the large decrease in area 
i-Q ' cannot be explained in terms of modification of the blackbody spectrum by electron scattering in the atmosphere 

of the neutron star. The height of the emitting region on the non-flaring neutron star is shown to agree with 
the height of the inner radiatively-supported accretion disk as found for sources in the ASCA survey of LMXB 
of Church & Balucihska- Church (2001). The decrease in height during flaring is discussed in terms of possible 
models, including radial accretion flow onto the stellar surface and the theory of accretion flow spreading on the 
neutron star surface of Inogamov & Sunyaev (1999). We demonstrate that the intensity of the broad iron line at 
6.4 keV is strongly correlated with the luminosity of the blackbody emission from the neutron star, and discuss 
the probable origin of this line in the ADC. Finally, possible reasons for non-detection of a reflection component 
in this source, and LMXB in general, are discussed. 
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1. Introduction 

X 1624-490 is one of the most unusual members of the 
class of dipping Low Mass X-ray Binary (LMXB) sources 
exhibiting periodic dips in X-ray intensity. It is generally 
accepted that dipping is due to absorption in the bulge in 
the outer accretion disk where the accretion flow from the 
companion impacts (White & Swank 1982). X 1624-490 
has the longest orbital period of the dipping sources at 
^21 hr. Dipping is deep, ~ 75% in the band 1-10 keV, and 
the source also exhibits strong flaring in which the X-ray 
flux can increase by 30% over timescales of a few thousand 
seconds (Church & Balucihska-Church 1995). Flaring has 
been known in LMXB for many years since it was recog- 
nized that ScoX-1 had quiescent states and active states 
in which strong X-ray flaring took place (Canizares et al. 
1975). White et al. (1985), using the Exosat GSPC, showed 
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that the continuum between 2-25 keV could be modelled 
by a blackbody plus unsaturated Comptonization spec- 
trum and that the flaring is caused by increases in the 
blackbody luminosity from 10% to 40% of the total. The 
blackbody temperature increased significantly in flaring 
together with some increase in emitting area. Hasinger et 
al. (1989) realized that ScoX-1 was a Z-track source, the 
source having different spectral states forming a Z-track 
on a hardness-intensity diagram with associated changes 
in QPO behaviour. It was also realized that several other 
bright sources exhibited similar flaring to that in ScoX- 
1. However, since then, little work has been carried out 
on flaring in LMXB, physical understanding of the flaring 
phenomenon is limited, and the transition between burst- 
ing in less bright LMXB and flaring in brighter sources is 
not understood. 

Investigation of X-ray dipping has however, been 
more extensive. The depth, duration and spectral evolu- 
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tion in dipping vary considerably from source to source. 
Moreover, spectral analysis of dipping sources provides in- 
formation not available in non-dipping sources. For ex- 
ample, spectral models are more strongly constrained by 
having to fit non-dip and dip data, thus showing clearly 
the nature of the emission, ft has been shown that the 
dipping LMXB sources are well described by assuming 
point-like blackbody emission from the neutron star plus 
Comptonizcd emission from an extended Accretion Disk 
Corona (ADC) (e.g. Church & Balucinska-Church 1995; 
Church 2001). In addition, the dip ingress and egress times 
can be used to obtain the sizes of the extended emission 
region if the absorber has larger angular extent that the 
source regions, as proven by dipping being 100% deep at 
any energy, which is often the case. This provides a ra- 
dius for the ADC, ^adc, typically 50,000 km or 15% of 
the accretion disk radius, although in X 1624-490, tadc is 
5.3 x 10 10 cm or 50% of the disk radius (Church 2001). 
There is also evidence that the ADC is "thin", having 
height-to-radius ratio of -10% (Smale ct al. 2001). This 
shows that the ADC is a very extended, flat region above 
the accretion disk, contrary to theoretical suggestions of- 
ten made that the Comptonizing region may be localized 
to the neighbourhood of the neutron star. 

Recent results of an ASCA and BeppoSAX survey of 
LMXB (Church & Balucihska- Church 2001) reveal a sur- 
prising result for the area of the neutron star emitting 
blackbody radiation for sources of different luminosities. 
When it is assumed that the emission region is an equa- 
torial strip of half-height h, it is found that there is good 
agreement between h and H, the half-height of the inner 
disk, assumed to be radiatively supported. It is not clear at 
the moment whether there is a direct causal link between 
these parameters, for example, if material flowed radially 
between the inner disk and the star, or whether there is an 
indirect link such that H is a measure of the total source 
luminosity Lxot and some physical process depending on 
Lxot determines the extent of the emitting area. In partic- 
ular, Inogamov & Sunyaev (1999) have shown that accre- 
tion flow regarded as a fluid meeting the neutron star at 
the equator is expected to spread vertically up the star to 
a height dependent on Lxot- A preliminary comaparison 
of this model with the survey results (Church et al. 2001) 
shows that observed blackbody luminosities are ~3 times 
larger than implied by the theory, suggesting that radial 
flow may be the dominant factor determining the black- 
body luminosity. The results of the present work will be 
viewed (Sects. 3 and 4) in terms of the h = H relation. 

X 1624-490 was previously observed in a long Exosat 
observation of 220 ks and with Ginga (Jones & Watson 
1989). The Exosat ME data revealed an apparently sta- 
ble lower level of dipping supporting the presence of two 
emission components, one of which was totally absorbed in 
deep dipping. Church & Balucihska-Church (1995) showed 
that the light curve at higher energies (> 5 keV) was 
dominated by flaring which can strongly modify the spec- 
trum. Thus if flaring occurs simultaneously with dipping, 
the investigation of spectral evolution in dipping will be 



difficult. By selecting sections of non-dip and dip data 
without apparent flaring, the above two-component model 
was found to fit the data well showing that in deep dip- 
ping the blackbody component was totally absorbed, and 
that the Comptonizcd component was relatively little ab- 
sorbed. However, with only one non-dip and one deep dip 
spectrum it was not possible to determine the extent of 
absorption of the Comptonizcd component. The Galactic 
column density of X 1624-490 is very high (~ 8x 10 22 
atom cm~ 2 ) so that a dust scattered halo of the source 
is expected, and Angelini et al. (1997) demonstrated an 
excess of the surface brightness above the point spread 
function in ASCA CIS. 

The 60 ks observation with BeppoSAX of 1999, August 
11 was almost free of flaring and fitting the radial in- 
tensity profile of the MECS instruments in four energy 
bands between 2.5 and 6.5 keV allowed a value of the op- 
tical depth to dust scattering of 2.4±0.4 to be determined 
(Balucihska-Church et al. 2000). Moreover, the very broad 
band of BeppoSAX allowed the Comptonization cut-off en- 
ergy to be determined for the first time, this having the 
value of 12+g 4 keV. Spectral evolution in dipping could 
be well-described in terms of an emission model consist- 
ing of point-like blackbody emission from the neutron star 
plus extended Comptonized emission from the ADC. It 
was shown that using a disk blackbody model instead of 
a simple blackbody resulted in an unphysical inner radius 
of 0.4 km, so this possibility could be discounted. Similar 
unphysical results were obtained for this model for the 
majority of the LMXB in the ASCA survey allowing us to 
associate the blackbody emission with the neutron star in 
LMXB in general. Additionally, the extended size of the 
ADC and high optical depth imply that all of the thermal 
emission from the disk will be reprocessed in the ADC. 
Thus, we do not expect any contribution to the observed 
thermal emission from the disk. 

Several RXTE observations have been made by the au- 
thors including observations on 1997, January 6-7, 1997, 
May 4, 6, 11 which failed to capture dipping. A 4-day ob- 
servation with RXTE was made in 1999, September span- 
ning 4 orbital cycles with deep dipping and extensive flar- 
ing. Results of analysis of the spectral evolution in dipping 
are presented in Paper I (Smale et al. 2001). The orbital 
period has not been refined since the Exosat determina- 
tion of 21±2 hr (Watson et al. 1985). In Paper I, an orbital 
ephemeris and improved period of 20.8778±0.0003 hr are 
presented based on 4.5 yr of RXTE ASM data. In addi- 
tion, a search failed to detect QPO at frequencies from 
Hz to kHz, which may be related to the high luminos- 
ity of the source. The much higher count rate and obser- 
vation duration in the 1999, September observation with 
RXTE revealed several new features including the domi- 
nance of electron scattering in the early stages of dipping. 
The main deep dips were surrounded by shallow "shoul- 
ders of dipping" in which a vertical shift of the spectrum 
relative to the non-dip spectrum was seen between 2.5 
and 25 keV. This demonstrated that electron scattering 
was dominant and photoelectric absorption was not de- 
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tectable indicating the highly ionized state of the outer 
absorber produced by the high luminosity of X 1624-490. 
A broad iron line at ~6.4 keV was first detected by Asai et 
al. (2000) using ASCA. Spectral evolution in dipping could 
be well described using the above two-component emission 
model having the basic form: AG*(ab* bb + pcf* (cpl + 
GAu)), where AG represents Galactic absorption, AB addi- 
tional absorption of the point-like blackbody in dipping, 
PCF the progressive covering of the extended Comptonizcd 
emission (described as a cut-off power law) and GAU the 
Gaussian line. Good fits were obtained by including the 
6.4 keV line in the same covering term as the Comptonizcd 
emission, showing that the line probably originates in the 
ADC. However, to the basic form were added complicat- 
ing terms for the effects of both electron scattering and 
dust scattering (Smale et al. 2001). 

In the present paper (Paper II), we present a study of 
flaring in X 1624-490 based on the observations of 1997, 
January and 1999, September during both of which strong 
flaring took place. Although flaring in ScoX-1 and similar 
sources has been known for many years, in which clear in- 
dications of increasing blackbody temperature were found, 
there has been no recent observational work aimed at un- 
derstanding the physical changes taking place, or the dif- 
ferences between bursting and flaring sources. In addition, 
flaring is clearly related to the flaring branch of the Z-track 
class of LMXB. X 1624-490 offers a unique opportunity for 
investigating these problems as it is both a dipping and a 
flaring source. 
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2. Observations 

X 1624-490 was observed using RXTE (Bradt, Rothschild, 

6 Swank, 1993) from 1999, September 27 18:22 UT - 
October 1 13:11 UT, for a total on-source exposure of 212 
ks. In addition to results from this main observation span- 
ning 4 orbital cycles of 20.87 hr, we present results from 
our observation of 1997, January 6 11:14 UT - January 

7 3:28 UT. The new ephemeris (Paper I) shows that this 
observation, intended to capture dipping, spanned phases 
0.11-0.89 where phase is dip centre, exactly missing deep 
dipping. However, there was a period of strong flaring and 
both observations are of interest because of the change in 
source luminosity between the observations, the 1997 ob- 
servation having a luminosity (1-30 keV) of 1.06 x 10 38 erg 
s" 1 , compared with 1.35 x 10 38 erg s _1 in the 1999 ob- 
servation, approaching closer to the Eddington limit for a 
1.4 Mq object. It will be seen that this difference appar- 
ently leads to significant differences in flare development 
(Sect. 3). For the analysis presented here of relatively slow 
flare development and decay, data from the PCA instru- 
ment in the Standard 2 mode were used having a time 
resolution of 16 sec. The PCA consists of five Xc pro- 
portional counter units (PCUs) numbered to 4, with a 
combined total effective area of about 6500 cm 2 (Jahoda 
et al. 1996). During the September 1999 observation only 
PCUs and 2 were reliably on throughout the observa- 
tion, whereas in 1997 all PCUs could be used. To allow 
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Fig. 1. PCA light curve of X 1624-490 with 32s binning, 
upper curve: the 1997, January observation in the band 
1.9-25 keV; lower curve: the 1999, September observation 
in the band 2.0-25 keV 

direct comparison of the lightcurves and spectra, results 
are presented here for both observations using PCUs 
and 2, the high count rate permitting this. Data were also 
obtained with the HEXTE phoswich detectors, which are 
sensitive over the energy range 15-250 keV (Rothschild 
et al. 1998). However, because of the relatively low break 
energy of the Comptonized spectrum at 12 keV (Paper 
I) producing a relatively low count rate in HEXTE, and 
because spectral evolution during flaring takes place pri- 
marily in the band 1-20 keV, the HEXTE data provided 
little useful additional information. 

Data analysis was carried out using the RXTE stan- 
dard analysis software, FTOOLS 5.0.1. Background sub- 
traction of the PCA data was performed using the 
"skyvle/skyactiv 20000131" models generated by the 
RXTE PCA team. The quality of the background subtrac- 
tion was checked firstly by comparing the source and back- 
ground spectra and light curves at high energies (50-100 
keV) where the source no longer contributes detectable 
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Fig. 2. PCA light curves of the 1999 observation of 
X 1624-490 in two energy bands: 1.7-5.0 and 10.0-20.0 
with 32 s binning 

events; and also using the same models to background- 
subtract the data obtained during slews to and from the 
source. We conclude that our background subtractions in 
the 2-20 keV energy range are accurate to a fraction of 
a count per second. The quality of both the background 
subtraction and the response matrices used was checked 
by including a careful analysis of Crab calibration obser- 
vations on 1999, September 26, October 13, November 8, 
and November 23, before and after our 1999 observation 
of X 1624-490 (see Paper I). Systematic errors of 1% were 
added to the spectra. 

3. Results 

3.1. The X-ray lightcurves 

Figure la and lb show light curves of the two observa- 
tions in the total PCA band. As the 1997 observation was 
made during PCA gain epoch 3, and the 1999 observation 
during epoch 4, we used appropriate PCA channel ranges 
to give very similar energy ranges in the two light curves: 
1.9-25 keV and 2.0-25 keV for 1997 and 1999, respec- 
tively. Figure 2 shows the 1999 observation in two energy 
bands: a low energy band (1.7-5.0 keV) in which dipping 
is predominant and a high energy band (10.0-20.0 keV) 
in which flaring dominates. 

In the 1997 lightcurve, a period of strong flaring is 
followed almost certainly by a shoulder of dipping, i.e. 
the period of shallow dipping that often precedes and fol- 
lows the main deep dipping. In Paper I, we showed that 
this consists of a reduction in intensity of the extended 
Comptonized emission by electron scattering in the outer 
layers of the absorber before the neutron star is overlapped 
by absorber and deep dipping commences. During the 
shoulder, the blackbody emission component was shown to 
have zero column density additional to the Galactic con- 
tribution; however, it is important to realize that this does 
not represent the quiescent, i.e. non-flaring, non-dipping 



emission which has a somewhat higher intensity. The 1999 
lightcurves in Fig. 2 demonstrate vividly that while dip- 
ping is essentially a low-energy phenomenon, flaring dom- 
inates the light curves at high energies, and can only, in 
fact, be detected above ~ 6 keV. The reason for this can 
be understood in terms of the non-flare blackbody tem- 
perature fcTeB ~ 1-3 keV which increases to ~ 2.2 keV in 
flaring, so that the spectra cross over at about 6 keV, re- 
sulting in a deficit below this energy and an excess above, 
compared with the non-flare spectrum. Close inspection 
of the 1999 light curve indicates that each flare is highly 
structured consisting of many shorter events, and we pro- 
pose here that flaring consists of a superposition of X-ray 
bursts occurring with much reduced separation compared 
with burst sources, due to the mass accretion rate being 
more than 10 times higher (at a luminosity of ~ 10 38 erg 
s- 1 ). 

3.2. Spectral evolution in flaring 
3.2.1. Selection of flare spectra 

Flare evolution is firstly displayed on X-ray hardness- 
intensity diagrams for both the 1997 and 1999 observa- 
tions in Fig. 3 (left panels). Hardness is defined as the 
ratio of intensities in the bands 9.8-24.9 keV and 4.9-9.8 
keV, and intensity is in the total band 1.9-24.9 keV. In 
the 1997 data (upper panel), strong flaring but only shal- 
low dipping took place, whereas in the 1999 data (lower 
panel), there was both strong flaring and deep dipping. 
In this case, dipping begins with intensity decreases of 
<~15% in an energy-independent way (horizontal part of 
Fig. 3) corresponding to a shoulder of dipping, and con- 
sistent with electron scattering of the non-dip spectrum 
in highly ionized outer layers of the absorber (Paper I). 
In deeper dipping, spectral evolution is complex, showing 
first a strong hardening as Comptonized emission is ab- 
sorbed, then a softening in deepest dipping (at < 80 c 
s _1 ) as the blackbody emission becomes totally removed 
leaving the residual spectrum harder. 

During flare evolution, systematic changes clearly take 
place in hardness-intensity very similar to the flaring 
branch of Z-track sources displayed on similar diagrams. 
Flare spectra selected for spectral analysis are also shown 
as average points for each spectrum in Fig. 3 (right pan- 
els). It was not initially clear that simply selecting flare 
spectra on the basis of intensity would be adequate to 
allow understanding of spectral evolution during flaring. 
Preliminary tests with this and other sources displaying 
track movement revealed that selection on intensity alone 
may produce scatter in the positions of the selected spec- 
tra about a line drawn through the centre of the track on 
the hardness- intensity plot. Such scatter leads to changes 
in spectral parameters not due to the development of flar- 
ing particularly in the Comptonized emission, and this 
complicates interpretation. Consequently it was decided 
to select spectra in such a way that the points lay along 
the centre of the track on the hardness-intensity diagram 
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Fig. 3. Flare evolution displayed on hardness-intensity diagrams; upper panel left: the 1997 observation; lower panel 
left: the 1999 observation. The dots show data that was excluded during selection of flare data, i.e. data in dipping, 
the shoulders of dipping, or flaring during the shoulders; open circles shows data accepted for flare analysis. The right 
hand panels in each case show the corresponding plots for spectra selected for analysis, allowing these to be compared 
with the total data (see text) 



that represents flare evolution. This was done for the 1997 
data by selecting data on the basis of both hardness ratio 
and intensity, defining essentially a box for each spectrum 
within Fig. 3. It can be seen that the spectra produced 
well represent the general track movement. In the case of 
the 1999 observation, spectra were first selected on the ba- 
sis of intensity alone, and when these data were mapped 
onto the hardness-intensity diagram, it was found that 
these points moved systematically along the flare evolu- 
tion track, obviating the need for selection on both hard- 
ness and intensity. This was a result of the much improved 
total count in spectra of the long observation. Comparison 
of results from this 'track-centre' technique with initial 
tests showed that the behaviour during flare development 
becomes systematic. 

For the 1997 observation, 9 flare spectra were se- 
lected using good-time intervals obtained from a hardness- 



intensity diagram made in the same energy bands as in 
Fig. 3 but using all 5 PCUs to reduce the errors on ev- 
ery point. Only the first 3 flares were used because the 
4th flare is in the shoulder of dipping. For the 1999 ob- 
servation, 12 flare spectra were selected in the same way 
avoiding dipping and shoulder-of-dipping data. The re- 
sults of fitting spectra selected in this way are discussed 
in the following section. Finally, we investigated whether it 
is adequate to add data from different flares by the proce- 
dure of intensity selection, and analysed data from single 
flares including the strong flare at ~180 ks in the 1999 
data. (Fig. lb). Results showed that this flare followed 
the same track on the hardness-intensity diagram as the 
total 1999 data, and spectral analysis results were also 
very similar. Consequently there was no need to restrict 
analysis to individual datasets. 



6 



Baluciriska-Church ot al.: Flaring in X 1624-490 



3.2.2. Spectral evolution during flaring 



Table 1. Ability of various spectral models to fit the 1999 
flare spectra simultaneously. The models are variations on 
the two-component bb + CPL model, except for the first, 
in which a second blackbody is added to attempt to model 
flaring 



Model 


X 2 /dof 


Additional blackbody 


1038/589 


Constant cut-off power law 


813/568 


Constant blackbody 


1493/568 


Constant blackbody radius 


904/557 


Constant blackbody temperature 


1283/557 


Constant power law index 


445/557 


Constant power law normalization 


370/557 


Free fitting 


354/546 



In Paper I, it was shown that non-dip, non-flaring emis- 
sion can be well-described by a two-component model con- 
sisting of point-like blackbody emission from the surface 
of the neutron star plus extended Comptonized emission 
from the accretion disk corona represented by a cut-off 
power law, plus a line at 6.4 keV. Dust scattering is signif- 
icant in this source, since photons from non-dip emission 
are scattered into the beam in dipping due to the time- 
delay in scattering. As the amount lost from the beam 
depends on the dip intensity; there is an excess due to the 
difference of these terms which makes a (non-dominant) 
contribution to the residual count rate in deep dipping. 
Outside dipping, the two terms balance and no effect is 
expected, and the spectrum is well-modelled by the two 
continuum components above. The cross-section for dust 
scattering has a strong energy dependence ~ E~ 2 there- 
fore the limited effects that are seen are at low energies 
and flaring at energies above 6 keV should not be influ- 
enced by dust scattering. In Paper I we exhaustively tested 
a range of one-component and two-component models us- 
ing non-dip and several dip spectra, which showed that 
only the above model was capable of fitting all of these 
spectra. Consequently in the present work we use that 
model to investigate spectral evolution during flaring, i.e. 
the model AG*(bb + CPL + GAu) where AG represents 
Galactic absorption, bb is the neutron star blackbody, CPL 
is the cut-off power law representation of Comptonization 
in an extended ADC and GAU is the Gaussian line. 

This model was firstly used to test various possibilities 
for flaring evolution, such as the Comptonized emission 
remaining constant, by fitting the 12 spectra from 1999 si- 
multaneously, and results are shown in Table 1 . The prob- 
ability of any model being correct was vanishingly small 
except for the last three, including attempting to model 
flaring by an additional blackbody component. The last 
model with all parameters free was more probable and ac- 



ceptable at >99% probability, and so this model was used 
in fitting the 1997 and 1999 observations. 

Spectral fitting was carried out by simultaneously fit- 
ting all flare spectra in each of the observations using 
the two-component continuum model described above to- 
gether with a line at 6.4 keV with width a fixed at 0.4 keV 
(as in Paper I) to prevent the tendency of a broad line to 
absorb neighbouring continuum. Initially, the Galactic col- 
umn density was free, however, this is not well-determined 
in RXTE and with free-fitting, values tended to be- 
come substantially higher (50%) than the well-determined 
value from BeppoSAX of 8.6±1.0xl0 22 atom cm~ 2 
(Balucihska-Church et al. 2000). Consequently, this value 
was fixed in the fitting. Using BeppoSAX with a spectrum 
extending to higher than 100 keV, it was also possible to 
obtain a good measurement of the Comptonization cut-off 
energy, while in the RXTE PCA with more restricted en- 
ergy band it is difficult to get reliable values, particularly 
of possible changes in this parameter during flaring. We 
show below that the most dramatic changes in flaring take 
place in the blackbody component; nevertheless, we tested 
the effect of different approaches to the cut-off energy as 
follows. Firstly, we allowed this parameter to be entirely 
free, or free but chained between individual flare spectra in 
each observation. Secondly, we fixed the value for all spec- 
tra at the best-fit value obtained from BeppoSAX of 11.8 
keV. Finally, we fixed the value at various values between 
11.8 keV and 100 keV. The quality of fit deteriorated as 
the value was increased, and maximum values of Eco of 
12 keV (1997) and 21 keV (1999) were found, after which 
X 2 /dof became unacceptable. With the cut-off energy free 
but chained between spectra, values of ~5 keV (1997) and 
^9 keV (1999) were obtained. However, the results pre- 
sented below were found not to be sensitive to Eco> par- 
ticularly the blackbody temperature and radius, and the 
systematic changes in these during flare development; con- 
sequently it was decided to adopt the approach of fixing 
the cut-off energy at the BeppoSAX value. Spectral fitting 
results are shown in Table 2 for the 1997 observation and 
in Table 3 for the 1999 observation. In Fig. 4 we show the 
blackbody radius Rbb, the temperature /cTbb, luminosity 
Lbb and the power law photon index T as a function of 
Lxot, for the 1997 and 1999 observations, respectively. 

In both observations, it was clear that strong changes 
take place during flaring in the blackbody component 
with A;Tbb increasing from ~ 1.3 keV to more than 2 
keV. Remarkably, the blackbody radius (where the emit- 
ting area is 47ri? BB ) decreased strongly in both obser- 
vations. In the 1997 data, Lbb remained approximately 
constant during flaring, and similar behaviour was seen 
in the 1999 data for L To t < 1-5 x 10 38 erg s -1 (equal to 
the largest Lxot in 1997); at higher Lxot, ^bb increased 
to ~ 3.5 x 10 37 erg s _1 at the peak of flaring. Thus, it is 
clear that over much of flare development, the decrease of 
Rbb offsets the increase of &Tbb 80 that Lbb remains con- 
stant. A consequence of the constancy of Lbb is that the 
increase in total luminosity is entirely in the Comptonized 



Baluciriska-Church ct al.: Flaring in X 1624-490 



7 



Table 2. Spectral fitting results for the nine 1997 RXTE 
PCA spectra fitted simultaneously; x 2 /dof = 400/405. 
The first row shows the quiescent source, and subsequent 
rows show the development of flaring. 90% confidence er- 
rors are given: luminosities are in units of 10 37 erg s -1 



Rossi-XTE PCA 



X 1624-490 



kT BB 


Rbb 


Lbb 


r 


LcPL 


keV 


km 


org s -1 




org s -1 


1.34 ±0.02 


10.2 ±0.4 


4.3 ±0.2 


2.09 ±0.06 


7.6 


1.40 ±0.03 


8.9 ±0.4 


4.0 ±0.2 


2.04 ±0.06 


8.1 


1.45 ±0.03 


8.3 ±0.4 


4.0 ±0.2 


1.98 ±0.06 


8.3 


1.54 ±0.05 


7.3 ±0.5 


3.8 ±0.3 


1.90 ±0.07 


8.8 


1.59 ±0.06 


6.7 ±0.6 


3.8 ±0.3 


1.88 ±0.08 


9.6 


1.78 ±0.08 


5.4 ±0.5 


3.8 ±0.3 


1.91 ±0.10 


9.8 


1.80 ±0.07 


5.4 ±0.5 


3.8 ±0.3 


1.85 ±0.08 


10.0 


1.90 ±0.07 


4.8 ±0.4 


3.9 ±0.3 


1.78 ±0.07 


10.9 


1.98 ±0.10 


4.5 ±0.5 


4.1 ±0.4 


1.79 ±0.10 


11.1 



emission. The power law photon index T decreased sys- 
tematically by 15% during flaring as seen in Fig. 4. 

The iron line had a non-flaring EW of 60± 20 eV in 
the 1997 data, and 90± 30 eV in the 1999 data. We do 
not show EW in Tables 2 and 3 for the growth of flares, 
because the continuum flux is changing; the line intensity 
is more relevant. In the 1997 data, the intensity of the 
6.4 keV line did not change, and it will be relevant that 
Lbb was also constant. In the 1999 data, the line inten- 
sity Juno did change and is plotted against Lbb and Lcpl 
(the luminosity of the Comptonized emission) in Fig. 5. 
This shows that during the initial growth of flaring the 
line intensity does not change significantly although Lcpl 
is increasing; Lbb is constant at this stage. However, in 
the second stage of flare growth, L B b increases (Fig. 4), 
and then the line intensity increases. In Fig. 5, the best-fit 
line to the variation of Li ne with L B b is shown, this tak- 



Table 3. Equivalent results for the 1999 data. x 2 /dof = 
354/546 for 12 spectra fitted simultaneously; first row is 
non-flare as before 



kT BB 


Rbb 


Lbb 


r 


Lcpl 


keV 


km 


erg s _1 




erg s _1 


1.37 ±0.04 


7.1 ±0.5 


2.3 ±0.2 


2.20 ±0.04 


11.2 


1.37 ±0.03 


7.1 ±0.4 


2.3 ±0.2 


2.17 ±0.03 


11.4 


1.42 ±0.03 


6.7 ±0.4 


2.4 ±0.2 


2.13 ±0.03 


11.5 


1.50 ±0.03 


5.9 ±0.3 


2.3 ±0.1 


2.12 ±0.03 


11.9 


1.62 ±0.04 


5.1 ±0.3 


2.3 ±0.1 


2.09 ±0.03 


12.1 


1.71 ±0.05 


4.5 ±0.3 


2.2 ±0.1 


2.05 ±0.03 


12.3 


1.89 ±0.06 


3.8 ±0.3 


2.3 ±0.1 


2.04 ±0.04 


12.6 


2.02 ±0.06 


3.4 ±0.2 


2.6 ±0.2 


2.01 ±0.04 


12.5 


2.16 ±0.06 


3.1 ±0.2 


2.8 ±0.2 


2.01 ±0.15 


12.7 


2.18 ±0.06 


3.1 ±0.2 


2.8 ±0.2 


1.95 ±0.05 


12.8 


2.22 ±0.07 


3.1 ±0.2 


3.0 ±0.3 


1.94 ±0.06 


12.8 


2.25 ±0.10 


3.2 ±0.4 


3.4 ±0.3 


1.96 ±0.11 


13.0 
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Fig. 4. Upper: Spectral fitting results from the 1997 ob- 
servation: blackbody parameters and power law photon 
index T. The blackbody luminosity is in units of 10 36 erg 
s -1 . Lower: Corresponding results from the 1999 observa- 
tion 



ing into account the 90% confidence errors in both y and 
x. This line has gradient 8.57 ± 0.69 x 10~ 5 photon cm~ 2 
erg -1 , showing that the line intensity is correlated with 
the blackbody luminosity with high significance. This is 
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Fig. 5. Variation of iron line intensity (in photon cm~ 2 
s _1 ) with the luminosity of the Comptonized emission 
(upper panel), and with the blackbody luminosity (lower 
panel). The best-fit line is shown in the lower panel (see 
text) 



the first time that a line variation has been linked to the 
variation of a particular continuum emission component; 
this will be discussed further in Sect. 4. 

The data of Fig. 4 are replotted in Fig. 6 as fcTeB 
against Lxot for both RXTE observations, plus a point 
from the BeppoSAX observation (Balucinska-Church et al. 
2000), when the source was substantially less bright. We 
have made comparisons of RXTE PCA Crab data in the 
band 3-20 keV (obtained before and after the 1999 ob- 
servation of X 1624-490) with BeppoSAX MECS data of 
1998, October, which shows that systematic differences be- 
tween flux from the two instruments amount to no more 
than 17%, (higher in RXTE). Fig. 6 shows that at the 
peak of flaring in the 1999 data, UTbb saturates at a con- 
stant value at the same time that Rbb becomes constant 
at a low value (Fig. 4). During this stage of flaring, Lbb 
increases, eventually becoming fixed when the emitting 
area and UTbb become constant. In Fig. 7, Lbb is plotted 
against blackbody temperature UTbb- The approximate 
constancy of Lbb m the early stages of flaring is appar- 
ent. Also shown is a line drawn through the penultimate 



point of the 1999 observation before the peak of flaring, 
representing the relation Lbb = AaT 4 . In this second 
stage of flaring as i?BB becomes constant, further increase 
in Lbb is due to temperature increase at ~ constant area, 
so that the data obey Stefan's Law. 



Rossi-XTE PCA 



X 1624-490 



XTE Jan 97 XTE Sept 



SAX Aug 98 



100 

, (1 - 30 keV; 10 3 



150 



Fig. 6. Variation of blackbody temperature UT-qb with 
source luminosity for the 1997 and 1999 RXTE observa- 
tions and the 1998 BeppoSAX observation 



Rossi-XTE PCA 



X 1624-490 




Fig. 7. Blackbody luminosity as a function of blackbody 
temperature which increases as flares evolve. The con- 
stancy of Lbb over most of the flaring can be seen, and 
the final increase of luminosity agreeing with a T 4 law, as 
the blackbody emitting area becomes constant (see text) 

Next, the results are re-expressed in the way used by 
Church & Balucihska-Church (2001) in a survey of LMXB 
sources using ASCA and BeppoSAX. In Fig. 8 we show 
the variation of the half- height h of the emission region on 
the neutron star as a function of Lxot- It is assumed that 
emission is from an equatorial strip on the star with area 
AnhR (a sphere intersected by two parallel planes) equal 
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Fig. 8. Evolution of the blackbody emitter during flaring 
in X 1624-490; the curve shows the equilibrium height of 
the radiatively-supported inner disk H eq 

to 4ttRq B . Thus, h = i?| B /i?, and it is assumed that R 
is 10 km. Also shown are calculated values of the equilib- 
rium half- height H eq of the inner radiatively-supported ac- 
cretion disk given approximately by the equation (Frank, 
King & Raine, 1992): 



H 



8 7rm p c 



r 



1/2- 



cq 



r 



1/2' 



(1) 



This equation shows that the disk height increases rapidly 
with radial distance from zero at the surface of the neu- 
tron star towards the equilibrium value H eq . Most of the 
increase takes place between r = 10 and r = 20 km, i.e. 
within 10 km of the stellar surface. In weak sources, with 
L~ 5 x 10 36 erg s -1 , radiative support is over a limited 
radial extent so that H eq is never achieved. However, in 
bright sources with L ~ 10 38 erg s _1 , such as X 1624-490, 
the equilibrium height is achieved, and the radiatively- 
supported disk extends to a radial distance of ^400 km. 
For the sources in the LMXB survey there was found to 
be an agreement h = H spanning 3 orders of magnitude 
in luminosity. 

In the case of X 1624-490, we calculated H eq from Lxot 
using the expression Lxot = GM M/R, where M is the 
mass of the neutron star taken to be 1.4Mq. However the 
increase of luminosity in flaring probably does not mean 
an increase in M (as in X-ray bursts). Thus the curve 
of H cq shown in Fig. 8 is appropriate for comparison of 
the non-flaring source in the 3 observations, but should 
not be taken to imply M increases in flares. It can be 
seen that there is good agreement between h and H eq for 
the BeppoSAX and 1997 RXTE observations. In the 1997 
data, h = 10 km, i.e. the whole neutron star is emitting. 
During flaring, there is a strong depression of h below 
H eq by a factor of ^5 at the peak of flaring. In the 1999 



observation, it can be seen that the non-flaring source has 
h already depressed below H eq . 

As flaring develops in the 1997 observation, Lxot in- 
creases from 1.15 x 10 38 erg s _1 to 1.5 x 10 38 erg s -1 . The 
1999 RXTE data shows similar behaviour, except that the 
non-flaring luminosity was higher at 1.3 x 10 38 erg s _1 
rising to 1.7 x 10 38 erg s _1 in flares. It is important to 
decide whether the initial depression of h below H cq is 
real so the effects of uncertainty in the Comptonization 
cut-off energy discussed previously were tested. The value 
of Eqo was allowed to vary between its smallest value 9 
keV obtained in free fitting and the maximum value 21 
keV. The effect of this was that the total luminosity (1-30 
keV) varied as might be expected, from the best-fit value 
of 1.3 x 10 38 erg s" 1 for E co = H-8 keV by ±0.1 x 10 38 
erg s _1 for the non-flaring spectrum. However, there was 
essentially no change in the height of the emitting region 
on the neutron star. We can conclude that the dramatic 
changes shown in Fig. 7 are not affected by this uncer- 
tainty. 



4. Discussion 

The main result of our analysis is a neutron star black- 
body radius which decreases strongly as flares develop. 
We firstly consider whether this can be influenced by the 
fact that we fitted a simple blackbody when it is often as- 
sumed that X-ray burst spectra, particularly at the peaks 
of bursts can have a modified blackbody spectral form, 
due to the effects of electron scattering in the neutron star 
atmosphere. There were two types of evidence for a mod- 
ification taking place: blackbody temperatures of 2.5-3.0 
keV implying super-Eddington luminosities in bursts, and 
an anti-correlation of blackbody radius Rbb with fcTeB 
in a particular source (below) when a simple blackbody 
was fitted. It should be mentioned at this stage that it 
would be unwise not to fit a simple blackbody model to 
the present data as the spectral form depends on one pa- 
rameter only fcTBB, allowing straightforward interpreta- 
tion of the results. Particular modified blackbody models 
have a flux depending on kT and electron density and so 
provide no information on the emitting area unless some 
assumption is made about the electron density. 

In the study of X-ray bursts, it was first shown by 
Swank et al. (1977) that the burst spectra were well de- 
scribed by a simple blackbody. This was supported by 
the fact that during burst decay, the flux varied approx- 
imately as T 4 (Hoffman ct al. 1977). By fitting a sim- 
ple blackbody model, burst peak blackbody temperatures 
of ^2 keV are generally found; however higher tempera- 
tures of ~3 keV have been measured, also obtained using 
high quality data from the Rossi-XTE PCA (Swank priv. 
comm.), so that high temperatures at the peak of some 
bursts cannot be doubted. As these high temperatures im- 
ply super-Eddington luminosities, theory was developed in 
which electron scattering dominated over free-free absorp- 
tion in the opacity of the neutron star atmosphere dur- 
ing bursts, and the surface value of the radiation source 



10 



Baluciriska-Church ct al.: Flaring in X 1624-490 



function departed from Planck form and had a modified 
blackbody form. Modified blackbody spectral forms have 
been given for a constant density atmosphere which we 
will designate mbbI (van Paradijs 1984) and for an ex- 
ponential atmosphere mbb2 (Zel'dovich & Shakura 1969). 
Modelling of the neutron star atmosphere has been car- 
ried out by Fujimoto et al. (1981) who derived conditions 
for stable and unstable hydrogen and helium burning, by 
Paczyhski (1983), by London et al. (1984, 1986), Ebisuzaki 
et al. 1984, and Ebisuzaki 1987. London et al. modelled 
the atmosphere numerically allowing for the effects of 
Comptonization. Although the simpler models mbbI and 
mbb2 have a spectral form which differs from the Planck 
form, and could, in principle, be detected if sufficient count 
statistics were available, the more sophisticated numerical 
models have forms which are similar to the Planck form 
although the peak energy has a different relation to kT. 
However, when close to the Eddington limit, a substantial 
low-energy hump is expected (Madej 1991). In our work 
on LMXB, we have never found any indication of a de- 
parture from simple blackbody spectral shape in the qui- 
escent spectra, all sources studied being well-fitted by us- 
ing a simple blackbody plus cut-off power law model (e.g. 
Church & Balucihska-Church 1995; Balucihska-Church et 
al. 1999; Smale et al. 2001). X-ray burst spectra have been 
well-fitted by a simple blackbody, e.g. using the Rossi- 
XTE PCA with its high sensitivity, it has been possible 
to examine burst spectra accumulated over times as short 
as 0.1 s, and these are well-fitted by a simple blackbody as 
in 4U 1728-34 (Swank 2001). One possibility remaining is 
that bursts with kT ~3 keV may have modified spectra, 
whereas those with kT ~2 keV do not. 



f-f absorption dominates 



0.5 



electron scattering dominates 



kT (keV) 



Fig. 9. Relative contributions of absorption and electron 
scattering to the opacity in the neutron star atmosphere 
as a function of number density and temperature kT 



Modification of the blackbody spectrum was expected 
on the grounds that in an X-ray burst, the density was 
"low" and the temperature "high", and we next address 
the question of how well the electron density is known. 



It is convenient to use a density-temperature diagram on 
which we can delineate the regions in which free-free and 
electron scattering processes dominate as shown in Fig. 
9. Density can be expressed as electron number density, 
as mass density p, or as the vertically integrated mass 
column density S. The problem is that the value of density 
used in various theoretical treatments of X-ray bursting 
are very different. Fujimoto et al. (1981) provide values 
of AM ~ 10~ 12 -10~ n M for the total mass needed to 
be accumulated on the surface of the star for a nuclear 
flash to occur. From this, £ = 10 8 -10 9 g cm~ 2 and if 
we assume the thickness of the atmosphere in the range 
10-100 cm, the average p is 10 6 -10 8 g cm~ 3 . However, 
the density can vary considerably from the mean, and be 
different in the hydrogen and helium burning shells (Joss 
& Rappaport 1984). London et al. (1984) assume a value 
of £ = 100 g cm -2 . Ebisuzaki (1984) gives a range of mass 
densities in the surface layer for a burst luminosity of 20% 
of Eddington from 1-10 4 g cm -3 , which correspond to £ in 
the range 5-5 x 10 4 g cm~ 2 . Paczyhski (1983) used values 
of surface mass (AM) similar to those of Fujimoto et al. 
(1981). Thus these values correspond to £ spanning the 
range 5-10 9 g cm~ 2 and to electron number densities n e in 
the range 10 24 -10 32 cm -3 . While an atmosphere with n e 
= 10 24 cm~ 3 is clearly dominated by electron scattering, 
this is not the case for higher values of electron density. 
The calculation of densities in neutron star atmospheres 
during bursts is complex and it does not appear possible 
yet to be definite about values. 

The second piece of evidence for modification of the 
blackbody in X-ray bursts was the apparent decrease in 
blackbody radius with increasing temperature kT when a 
simple blackbody model was fitted to Tenma and Exosat 
data of bursts in XB 1636-536 (Inoue et al. 1984; Sztajno 
et al. 1985). In these data, Rbb decreased by about 40% 
as fcTeB increased from 1 to 2 keV. Sztajno et al. showed 
that the apparent decrease was consistent with a constant 
emitting area if the model atmospheres of London et al. 
(1984, 1986) were used. We will make the same compari- 
son with the present data. London et al. construct mod- 
els for various values of the emitted flux /, which they 
re-express in terms of an effective temperature T c ff, de- 
fined via f ~ T 4 ff . Fitting a blackbody to the emission 
of their model atmospheres provides a spectral or colour 
temperature T spec , so that if the luminosities of the fitted 
blackbody and generated spectrum are equal: 



4 tt Rl B a T 4 _ = const ■ R 2 T% 



spec 



cfr 



(2) 



Assuming R is constant gives -Rbb °c (7efT/T spec ) 2 . Using 
the spectral hardening factor T S pcc/7cff ~ T^g 5 of London 
et al., we obtain i?BB ~ l/Tsp 7 ^ for the expected depen- 
dence of the apparent decrease in blackbody radius with 
increasing measured burst temperature. In the case of the 
present flare data if we assume some modification of the 
spectrum takes place, as kT increases from 1.34 to 1.98 
keV (1997 data), we would expect a decrease in Rbb of 
34%, whereas the actual increase is by a factor of 2.3, and 
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Fig. 10. Test of the hypothesis that the apparent decrease 
in blackbody area in flaring is not real by comparing 
the 1997 results with the R-T relation predicted by the 
London et al. model atmosphere (dashed line) assuming 
that the emitting area remains constant (see text) 



these changes are compared in Fig. 10. Thus, it does not 
seem very likely on the basis of the London et al. model 
that the emitting area can actually remain constant while 
appearing to decrease by this large factor. 

If we use the simpler mbbI modified blackbody model, 
the luminosity is related to the electron density via 
const ■ R 2 (fcT) 9 / 4 n\ (with a similar relation for mbb2). 
To demonstrate the constraints imposed on n c , we fitted 
the 1999 data using the model mbbI + CPL. This fit- 
ting provided good values of x 2 /dof; results are shown in 
Table 4. We assume that the spectrum of the non-flaring 
source is an unmodified blackbody, but that mid-flare and 
peak-flare can be modified. In the non-flare emission, the 
blackbody radius Rbb was 7.1 km. Assuming that this 
does not change in flaring, we use the measured mbbI 
normalization and kT values to derive the electron density 
at the peak of flaring, n e = 1.9 x 10 23 cm -3 . Not only is 
this fine-tuning of electron density values, but the value is 
relatively low in terms of expected values in model atmo- 
spheres as discussed above. Moreover, the electron density 
would have to decrease by a factor of two between mid- 
flaring and the peak of flaring to keep R constant. It is 
unlikely that n c decreases in flaring and so for this simple 
model also, it is unlikely that the emitting area does not 
decrease. 

However, it is not proven that a modification does take 
place in flaring, and there are arguments that it does not. 
Firstly, the value of kT ~2 keV obtained at the peaks 
of flares does not demand a modification of the spectrum 
in the way that super-Eddington values of kT = 3 keV 
do in X-ray bursts. Secondly, X-ray flaring is much less 
well understood than X-ray bursting and so our knowl- 
edge of densities in the neutron star atmosphere is even 
more uncertain as no theoretical model of the atmosphere 
during flaring exists at present. We have previously sug- 



gested (Sect. 3.1) that an X-ray flare consists of a su- 
perposition of X-ray bursts. The recurrence time between 
these is expected to be short as observed, because of the 
much higher values of mass accretion rate compared with 
a burst source. The 1-30 keV luminosity of X 1624-490 
in the observations discussed here is ~ 10 38 erg s _1 , i.e. 
considerably higher than that of a typical burst source. 
Consequently it is likely that the densities in the atmo- 
sphere are more than 10 times higher making it less likely 
that electron scattering will dominate the opacity. In sum- 
mary, there is no observational evidence of a modification 
of the blackbody spectrum. However, if a modification did 
take place, it is unlikely that this can explain the large 
observed decrease in area, and we conclude that the main 
results of the analysis made using a simple blackbody are 
correct. 




Non-flaring 




Fig. 11. Schematic diagram showing the result from spec- 
tral fitting that the emitting area of the neutron star 
blackbody decreases strongly during flaring. The top panel 
shows the calculated height H of the radiatively-supported 
inner disk which agrees well with the neutron star emitter 
height h in the non-flaring state. One possibility (see text) 
is that the vertical extent of the inner disk is decreased in 
flaring (lower panel) 

We now discuss possible reasons for the decrease in 
emitting area. Results of the ASCA survey of LMXB 
(Church & Balucinska-Church 2001) showed that the level 
of neutron star blackbody emission can be explained by 

Table 4. Results of fitting non-flare emission by simple 
blackbody and flaring by a modified blackbody spectral 
form, norm has units of 10 37 erg s _1 at 10 kpc; electron 





Non-flare 


Mid-flare 


Peak-flare 




Simple BB 


Modified BB 


Modified BB 


kT (keV) 


1.37±0.04 


2.30±0.07 


3.46±0.10 


norm 


1.02 ±0.09 


1.86 ±0.11 


3.35 ± 0.53 


R (km) 


7.1±0.5 


7.1 


7.1 


n c 




3.7 ±0.7 


1.9 ±0.3 
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assuming that the half-height of emitter on the star h 
is equal to the half-height H of the inner radiatively- 
supportcd disk (see Sect. 1). Two possible explanations 
exist (see Sect. 1): either there is radial flow across the 
gap between disk and star (mapping the value of H onto 
h), or H is simply a measure of M, and the mechanism 
of Inogamov & Sunyaev (1999) (accretion flow spreading 
on the stellar surface) determines h. If the first of these is 
correct, then radiation pressure can affect the inner disk 
by material being removed at the upper and lower sur- 
faces by radiation pressure, leading directly to a decrease 
in h and of M reaching the star in flaring, as shown in 
Fig. 11. If the second mechanism is correct, then we need 
the accretion flow between disk and star to be impeded 
by radiation pressure. In cither case, it is important to 
compare the levels of blackbody emission obtained from 
spectral fitting with the Eddington limit to test whether 
the effects of radiation pressure should be expected. If the 
emitting area is reduced at the peak of flaring to a narrow 
equatorial strip, radiation will not be isotropic. For an ob- 
server at angle 9 from the equatorial plane, the observed 
flux will be fobs — f ■ cos# where / is the flux in the plane. 
For X 1624-490, 9 -30°. To test whether Eddington limit- 
ing takes place in the equatorial plane, we should compare 
/ = /obs/cos0 with / Edd . 

Taking into account the non-Newtonian gravitation 
near the neutron star and possible variations in the com- 
position of neutron star atmosphere so that the opac- 
ity may not equal the Thompson scattering opacity, 
there will be a local value of critical luminosity L cr ;t = 
-C'Edd • ^ • (1 + z) where L E dd is the normal Eddington 
limiting luminosity, kt is the Thomson scattering opac- 
ity, k is the actual opacity and z the gravitational redshift 
(Paczyhski 1983). For a distant observer, this critical value 
will be reduced by the gravitational redshift to L C rit,oo, so 
that 

_ AirGMc K T 1 . , 

-^crit,oo — / -i . \ I'jJ 

KT K (1 + Z) 

where A-kG M c /kt is the standard Eddington limiting lu- 
minosity Lftdd — 1-76 x 10 38 erg s" 1 for a 1.4 Mq neutron 
star. This depends on composition since kt — 0.2(1 + X) 
where X is the fraction by weight of hydrogen. The factor 
1/(1 + z) = yjl - 2GM/Rc 2 in terms of the mass M and 
radius R of the star. We will compare our measured black- 
body luminosities with the Eddington limit by comparing 
with L C rit,oo, assuming simply hydrogen composition and 
a value of z = 0.34 appropriate to the surface at R = 10 
km of a 1.4 Mq neutron star. Thus £ C rit,oo 

= 1.31 x 10 3S 

erg s^ 1 . 

In both the 1997 and 1999 data, the observed black- 
body luminosity Lbb remains less than 4 x 10 37 erg s _1 . 
Thus the flux in the equatorial plane corresponds to a lu- 
minosity of ~ 5 x 10 37 erg s _1 , i.e. about a factor of 2.5 
less than L C rit,oc- This calculation suggests that the effects 
of radiation pressure are important. In the 1997 data, Lbb 
remains about constant at 4 x 10 37 erg s _1 and the black- 
body radius Rbb decreases from its non-flaring value of 



10.2 ± 0.4 km as fcTeB increases. In the 1999 data, Rbb 
starts at a somewhat lower level than in 1997 at 7.1 km 
so that the effects seen in flaring are already taking place 
in the non-flaring emission. The effects of radiation pres- 
sure clearly extend over a range of blackbody luminosity 
values, but at the peak of flaring further increase of Lbb 
is prevented due to both &Tbb an d -Rbb reaching stable 
values. Exactly how the lower limit of i?BB is caused is 
not clear. The processes discussed of increasing IzTbb and 
decreasing Rbb constitute a mechanism for limiting the 
blackbody luminosity and thus an Eddington limiting ef- 
fect not previously known. At the peak of flaring, the mass 
accretion flow to the star will be reduced below the non- 
flaring case by the changes in the inner disk and star, the 
radiation pressure will fall, and the height of the emitting 
region on the star increase. Thus, the flare is self-limiting 
and the system will return to the non-flare state on a sim- 
ilar timescale to that of flare growth. 

We next discuss further the implications of the rela- 
tion h = H obtained from the ASCA and BeppoSAX sur- 
vey (Church & Balucinska-Church 2001), and supported 
by the present work and the two possible explanations re- 
ferred to previously. The first explanation involves radial 
flow between inner disk and star. Although there has been 
a large amount of theoretical work on advective flow in 
black hole binaries, e.g. Abramowicz et al. (1996), and ref- 
erences therein, and on the nature of the disk itself, such as 
on slim disks (Abramowicz et al. 1988; Igumenschev et al. 
1998), there has been relatively little corresponding work 
on the accretion flow between the inner disk and neutron 
star in LMXB. It is not known whether the sonic point lies 
within the star or within the inner disk. However, Popham 
& Sunyaev (2001) investigate the nature of the inner ac- 
cretion disk around a neutron star, assuming a boundary 
layer within the disk, and show that the radial velocity 
increases by about two orders of magnitude leading to ad- 
vective flow within the disk and increase in disk height 
due to increased heating. This does not answer the ques- 
tion, of course, as to whether material can flow across the 
gap between the inner disk edge and the star, although 
Popham & Sunyaev comment that they would intuitively 
expect material to fall onto the star in a wide belt. Two- 
or three-dimensional modelling may reveal whether such 
radial flows are possible. It has been suggested that the 
model of radial flow between an inner radiative disk and 
the star is invalid as the inner disk may be unstable to 
the radiation pressure instability discussed by Janiuk et 
al. (2000) for a 10 Mq black hole. However, scaling from 
a 10 Mq black hole to a 1.4 Mq system indicates an in- 
stability timescale of ~30 s, and this, even if present, may 
not disrupt the disk. 

The alternative explanation of our results involves the 
flow of accretion material on the surface of the neutron 
star as proposed by Inogamov & Sunyaev (1999). The 
ageement between h and H would result from H being 
a measure of the total luminosity, not from H directly de- 
termining h in some way. However, preliminary results of 
a comparison of this model with the survey results show 
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that the observed blackbody luminosities are —3 times 
larger than predicted by the theory (Church et al. 2001). 
This implies therefore, that radial flow may be the domi- 
nant factor determining the neutron star blackbody emis- 
sion. The present results on flaring suggest for the radial 
flow explanation that radiation pressure can reduce the 
inner disk height so that both the emitting area on the 
star, and the mass accretion rate reaching the star are re- 
duced. Possibly, the accretion flow could be affected with- 
out the disk height changing. In the case of the Inogamov 
& Sunyaev mechanism, the decrease of M would lead to 
a decrease of emitting area. 

The present results are difficult to reconcile with the 
possibility that the blackbody emission in LMXB could 
originate in the accretion disk as multi-colour disk (MCD) 
blackbody emission. There would be no reason to expect 
any agreement between H and h, as h would be related to 
the inner radius n of the disk blackbody (via h = rf/R). 
The survey of LMXB (Church & Balucihska-Church 2001) 
showed that for the majority of sources, a Comptonization 
plus MCD blackbody model required unphysical values of 
inner disk radius, substantially less than the radius of the 
neutron star. The absence of disk blackbody emission in 
LMXB can however, be understood in terms of the very 
extended flat, nature of the ADC revealed by dip ingress 
time measurements (Sect. 1). An extended ADC of radius 
typically 50,000 km and large optical depth (Church 2001) 
will completely cover the hot X-ray emitting inner accre- 
tion disk, so that all thermal emission will be Comptonized 
and no disk blackbody component seen. 

Finally, we will discuss the broad iron line at —6.4 keV. 
In Paper I, we reported the detection of this line in non- 
dipping, non-flaring emission. It has an equivalent width 
of -60 eV and -90 eV for the 1997 and 1999 observations, 
respectively Asai et al. (2000) made the first detection of a 
broad line in this source at 6.41q; 3 keV using ASCA data; 
with an EW of 13+g 3 eV. In the present paper, we show 
that the line intensity is proportional to the luminosity of 
the blackbody component, the first demonstration that a 
line intensity is linked to a particular component of the 
continuum. Iron line emission is not well understood in 
LMXB with conflicting models for line origin. Asai et al. 
(2000) from their extensive study of iron lines in LMXB 
with ASCA find that the line energy had a mean value 
of 6.57 keV in 20 sources and argue that lines are formed 
by radiative recombination of photoionized plasma in the 
ADC, also proposed by Hirano et al. (1987) and White et 
al. (1985, 1986). In many sources, the line energy is —6.7 
keV strongly suggesting a blend of Fe XXV lines at 6.63, 
6.67 and 6.70 keV formed by recombination. Although 
there is some scatter in the line energies between extremes 
of 6.2 and 7.0 keV it is thought likely that the lines in all 
sources have a common origin. In general, measurement 
of line energies of 6.4-6.5 keV implying fluorescence has 
been taken to indicate origin in the accretion disk (e.g. 
Barret et al. 2000), although the number of detections of 
iron lines at these energies in LMXB is small. In the case 



of X 1624-490, the variation of the line in dipping (Paper 
I) is well described by giving the line the same covering 
factor as the ADC which supports the origin of the line in 
the ADC although the measured energy is unusually low 
at 6.4 keV. For the line energy to peak at 6.4 keV, an ion- 
ization parameter £ — L/nr 2 —100 is required which may 
be possible if the density in the ADC is high. However, 
it is difficult to reconcile this line with the theoretical ex- 
pectation that the electron temperature in the ADC of 
several keV might be expected to prevent formation of a 
6.4 keV line. 

An iron line produced in the ADC will have no bear- 
ing on the non-detection of a reflection component in this 
source. In fact, there has been a general lack of detec- 
tions of reflection components in LMXB and reported de- 
tections in a small number of sources have been ambigu- 
ous. Barret et al. (2000) presented evidence for a weak 
reflection component in GS 1826-238 and SLX 1735-269. 
However, in the first case, the spectrum was modelled with 
a one-component Comptonization term only but without 
a thermal component, which our previous work, e.g. the 
ASCA survey of LMXB (Church & Balucinska-Church 
2001) indicates is always present and would be detectable 
at the measured luminosity of the source. In the second 
case, the quality of fit was acceptable (x 2 /dof = 45/59) be- 
fore a reflection component was added. In a Ginga obser- 
vation of XB 1608-522, broad residuals above 7 keV could 
be modelled either by partial absorption or a reflection 
component (Yoshida et al. 1993). In all of our previous 
work using the bb + CPL model, there has been no trace 
of broad residuals indicating the presence of an additional 
continuum component. We can suggest a possible reason 
for the lack of detections. We have demonstrated the very 
extended, flat nature of the ADC (Church 2001; Smalc et 
al. 2001), with radius tadc typically 50,000 km, or 15% of 
the disk radius, and particularly large in X 1624-490 (see 
Sect. 1). Thus the ADC will cover all of the inner disk 
and act as a "hot" reflector preventing illumination of the 
disk by the neutron star so preventing reflection. In ad- 
dition, the large values of optical depth obtained for the 
ADC (Church 2001) will mean that any reflection compo- 
nent due to illumination by the ADC itself will not pene- 
trate the ADC to the observer. In CygX-1, tadc is much 
smaller (Church 2001) and there is evidence for a reflec- 
tion component, suggesting that LMXB neutron star bi- 
naries may differ from Black Hole binaries in this respect 
due to the strong effects of the neutron star in forming an 
extended ADC. 

Future work will concentrate on further investigation 
of flaring in sources such X 1624-490 and on the flaring 
branch of Z-track sources. It is clear that spectral evolu- 
tion during flaring can provide valuable information on the 
inner disk/neutron star interface not otherwise available. 
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